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In order to account for possible nonstatistical fluctuations in a hadronizing system while using a 
statistical approach, one has to resort to its nonextensive version. The new parameter q is intro- 
duced to be directly connected to the variance of an observable (with q = I one returns to usual 
statistics). We demonstrate how this approach allows composing fluctuations of different observ- 
ables and show that it leads to a specific sum rules proposing this to be verified experimentally. 
We discuss the ensembles in which all relevant quantities are allowed to fluctuate. By introduc- 
ing correlations between the observables, the relations connecting these variables are shown to 
generalize the so-called Lindhard thermodynamic uncertainty relations. This is illustrated using 
multihadron production data. We show that fluctuations from different components of collision 
phase space are correlated and that the strength of these correlations depend on the relevant q- 
parameter 
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1. Introduction 

Nowadays it is a standard procedure to use a statistical approach to model high energy mul- 
tiparticle production processes [|l|]. However, it has also been realized that data on many single 
particle distributions demonstrate a power-like behavior, rather than the expected simple exponen- 
tial one. In addition, multiparticle distributions are broader than naively expected. The proposition 
put forward some time ago is that, rather than invalidating a simple statistical approach, these 
observations call for its modification towards inclusion of a possible intrinsic, nonstatistical fluctu- 
ations, usually identified with fluctuations of the parameter T identified with the "temperature" of 
the hadronizing fireball (cf. reviews [|2|]). The introduced parameter q, known as the nonextensivity 
parameter, is shown to be directly connected to the variance of T , 

q = \ + (Oj with (Oj = ^j^2 ' (^-^^ 



and one obtains the so called Tsallis distribution 



, 2-^ / E\ 2-q 



^ 1 / E\ 
^-exp(--). (1.2) 



It reduces to the usual Boltzmann-Gibbs form for ^ = 1 (in which case one recovers the usual 
statistical model) ^ Such an approach is also known as Tsalhs statistics [^]. 

2. Results 

In this presentation we shall not go into the details of Tsallis statistics and its applications in 
the field of multiparticle production observed in high energy experiments (see ||^, |^). Here we 
concentrate on the problem of the coexistence and interconnection between fluctuations observed 
in different observables and, in particular, on the possible correlation between them which can, 
presumably, be measured experimentally. In any collision one observes a number of secondaries 
distributed in phase space {y,PT} (with rapidity y and transverse momentum pr)^. The multi- 
phcity distributions, P{N), rapidity distributions, dN/dy and distributions of transverse momenta. 



dN/dpT (or, transverse masses /I7 = y + p^) are usually measured. Among them, only P{N) 
refers to the entire phase space, the other are limited either to its longitudinal (y) or transverse (pj) 
parts. Therefore it is natural that nonextensivity parameters describing fluctuations present in these 
observables (visible either as a broadening of P(A'^) or power-Uke tails in other distributions) are 
different. This is clear from Fig. |l| where we show that, whereas fluctuations in full phase space 
depends strongly on energy, those confined to its transverse part depend on energy weakly. Fluc- 
tuations in the longitudinal part of phase space represent a special case as they strongly depend on 
data from which they are deduced and also on whether they are corrected for simultaneous fluc- 
tuations in pt or not [|l0|]. This is shown in Fig. |l| by the star symbols. As shown in [|lT]], all 

'This idea has originated in |^] and has been further developed in |^]. It forms a basis for so-called siiperstatistics 
|5|. The problems connected with the notion of temperature in such cases have been addressed in a recent book |^. Cf. 
1 2 1 for further references concerning applications of this approach. Recent examples of spectacular power-like behavior 
has been reported by PHENIX ^ and CMS |||| recently. 

^Where y = (1/2) In [{E + pi) / {E — pi)], with E being energy of particle and pi its longitudinal momentum. 
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these fluctuations can be connected when one considers an ensemble in which all variables char- 
acterizing an event, namely, the energy U, temperature T and multiplicity N, can fluctuate with 
relative variance co| = Var{X)/ < X >^, X = U, T, N. In this case, it occurs in a natural way 
that a correlation coefficient p 



Cov{U,T)/y/Var(U)Var{T) between the two of them (here U 
and T chosen) is necessary to fully describe the event. Generalizing the Lindhard thermodynamic 



uncertainty relations [12] one obtains [11 1: 
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Figure 1: Examples of energy dependence of the nonextensivity parameter q as obtained from different 
observables. Open symbols show q obtained from multiplicity distributions P{N) (fitted by ^ = 1 + 1 with 
l/k = — 0.104 + 0.029 Ini) [Q. Solid symbols show q = qT obtained from a different analysis of transverse 
momentum distributions, /{pr)- Data points are taken from, respectively, [|T^ for [Wibig], ]|^] for [CMS], 
from [ p^ ] for [NA49] (data on ixj) and from [ p6| ] for [UA5]. The dotted line represents a fit obtained in 
[0 (where qr = ^ -25 — O.SSi^*"*^^) and the full line comes from Eq. (p^ (for p = 0). Stars show q ^ qi 
obtained from dN/dy: open star taken from [|l6| (the same as used in [|17[]) shows the uncorrected q value, 
whereas the solid star indicates its corrected value [IlQ]. 



q-l 



(N) 



(0l, + (0j-2p(0u(0T- 



(2.1) 



What does the p tell us? For example, a large energy U (i.e., large inelasticity of reaction, K) can 
result either in a large number of secondaries of lower energies (for which p < 0) or a smaller num- 



ber of larger energies (for which p > 0, different models give different predictions). Eq. ( |2. Ij ) tells 
us that, in principle, the coefficient p is a function of all the nonextensivity parameters involved. 
Let us then specify the part of fluctuations of T which go in the transverse direction by a, and write 



qr 



1 



qi-\ = (ol + {\ 



(2.2) 
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From (2.2) one finds 

(oj = ^{qT-^), w^ = (^L-l)-(l-«)«| = (<?L-l)--^(<?r-l). (2.3) 



Then, from (2.3) and (2.1) one gets 



q-\ = {qL-\) + {qT-l)-2p(0uG>r ^ ^ - 1 = (^^ - 1) + (^r - 1) • (2.4) 



Assuming that 



Eq. (2.4) reads 



(0^ = K^(0^, (2.5) 



q-\ = {q^-\) + {gj--\)-2p^ {qr - 1) (2.6) 



whereas from (2.5) and (2.3) one has that 



That resuhs in a useful correlation coefficient given in terms of different (in principle measured) 
fluctuations: 

P= , ; a = ^. (2.8) 

i/afcUl 



Notice: 

• to get the correlation coefficient p, one has to know all the fluctuations, i.e., both in the entire 
phase space, q, as separately in its transverse, qr, and longitudinal, qi, parts; 

• out of them the best known is q (no corrections needed); 

• for qr the corrections are small and can be neglected; 

• however, for qi^ the corrections are large and must be accounted for (cf.. Fig. [l]). 

The example of feasibility of deducing p from data is presented in Fig. for data on p + p at 546 



GeV [16]. In this case from P(A'^) one gets q = 1.27, from distribution of pr one gets qr = 1.09 
whereas from the original qi = 1.48 one gets after correction qi = \A4 (cf. Fig. |]). 

To summarize, one can say that it is a priori possible to obtain p from experimental data, but - 
as shown here, this is difficult. Because of large estimation errors one cannot at this stage exclude 
the possibility that simply p = 0. Our preliminary result, cautiously indicating that, perhaps, p < 0, 
would prefer the production of a large number of particles of lower energies. However, statistically 
we cannot at this moment make a decisive statement. 



In [lllh we used a = 2/3 and )f = 1; for p = 0, the actual values of a and K parameters are irrelevant. 
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Figure 2: Left panel: the example of p obtained from Eq. (2.8). The shaded area shows the extent of the 
possible error, due to uncertainty in fixing qi- Right panel: the sensitivity of the correlation parameter p to 
a, specifying the fraction of fluctuations which goes into the transverse direction (in the case of full isotropy 
one has a ~ 2/3). 



3. Summary 

We discuss the hadronization process where all variables characterizing an event, namely the 
energy U, temperature T and multiplicity N fluctuate. In this case, it follows that correlation p 
between two of them (here chosen as U and T), is necessary to fully describe the event. The 
question is: what is its meaning and can it be estimated experimentally? The answer to the first 
is that this parameter tells us how the available energy is used: either for production of particles 
or, rather, for making them more energetic. In what concerns the second, our preliminary result 
(indicating that, perhaps, p < ) would prefer the first scenario, however, statistically we cannot 
say so at the moment. 



We conclude with two remarks. First, in what concerns Eq. (2.1), in the literature [18] there 
is similar relation comiecting the volume, V, pressure, P and temperature, T: (Op = (dy + (Oj [p^]. 
Secondly, when all variables, U , N and T fluctuate, the pairs of variables, {U,N) and {U, T), cannot 
aU be independent because 

Var{U) = {T)Cov{U,N) + {N)Cov{U, T) (3.1) 



(cf., [pi]]). This means that, in general, 

(Ou = p{V,N)(ON + p{VJ)(»T. (3.2) 



where p{X,Y) denotes the corresponding correlation coefficients between variables X and Y 
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'*For the completeness of presentation one must also mention that the notion of Tsallis statistics still remains a 
subject of hot debate (cf., for example, JT^]). Nevertheless, at the moment it is proved [^] that fluctuation phenomena 
can be incorporated into a traditional presentation of thermodynamics and result in many new admissible distributions 
satisfying thermodynamical consistency conditions among which is the Tsallis distribution (|l.2| ). We close this remark 
by noting that a recent generalization of classical thermodynamics to a nonextensive case pip clearly demonstrates its 
feasibility. 
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